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1 Introduction 



We consider the class of nonnegative weak solutions of the fractional diffusion equation 

dtu + {-Ay{u"') = in(0,T)xM^ (1.1) 
where m>0,0<,s<l,(i>l, and T > 0. The precise definition of the fractional Laplacian is given 



in Appendix 8.1 In most of the paper we assume that the initial data are given: 

u{0,x) = uo{x), (1.2) 

where in principle uq G L^(M'^) and uq > . However, in Section [2] we consider solutions of possibly 
changing sign and then we use the notation u"^ to mean for powers of signed functions. We refer 

to im [T2] for the basic theory of existence and uniqueness of weak solutions for the Cauchy problem 



(1.1) — (1.2). These papers also comment on the physical motivation and relevance of this nonlocal 
model, describe the main results on L'^ and C" regularity, and give references to related literature. 
Recently, the existence and properties of Barenblatt solutions for the Cauchy Problem was established 
in |18j . For s = 1 we recover the classical porous medium/fast diffusion equation, whose theory is 
well-known, cf. |17j . We will call the case s = 1 the standard diffusion case. 

The main purpose of the paper is obtaining quantitative a priori estimates of a local type for the solu- 
tions under consideration. Such estimates were obtained for the standard PME by Aronson-Caffarelli 
[U and by the authors for the standard FDE [41 El |6] . This is not always possible for the present model 
due to the nonlocal character of the diffusion operator, but then global estimates occur in weighted 
spaces. The results take different forms according to the value of the exponent m, a fact that is to 
be expected since it happens for standard diffusion. The case m > 1 is called the (fractional) porous 
medium case : contrary to the standard porous medium equation, it does not have the property of finite 
propagation, an important difference established in |ll[|12j . The range of exponents m G (0, 1) is called 
the (fractional) fast diffusion equation, and it has special properties when {d — 2s) /d =: mc < m < 1 , 
which we call the good fast diffusion range. When < m < mc it is known that some solutions ex- 
tinguish in finite time, which is a clear manifestation of the change of character of the equation, since 
solutions of the Cauchy problem exist globally in time and are positive everywhere in Q = (0, +00) x M'^ 
if 771 > m-c ■ 

Outline of the paper and main results. In Section [2| we derive upper bounds in form of weighted 
estimates, valid for nonnegative solutions of the Cauchy problem in the whole fast diffusion range 
< m < 1. Actually, they are valid for the difference of two ordered solutions, the precise statement is 



given in Theorem 2.2 Contrary to the purely local estimates known in the standard fast diffusion 
case, cf. [10], the estimates for s < 1 are valid in weighted L^-spaces and the weight must decay at 
infinity with a certain decay rate, not too fast, not too slow. This is again a manifestation of the 
nonlocal properties of the fractional Laplacian. The estimates will be important as a priori bounds for 
solutions, or families of solutions, through the rest of the paper. 

In Section [3] we use the estimates of Section [2] to construct solutions for initial data that belong to 
weighted L^-spaces, in particular for data < uq G L^^^CS.'^) such that uo{x) grows less than Odxp**/™) 
as |x| — )• 00, in particular for all bounded data. These solutions can be uniquely identified as minimal 
solutions in a precise sense and satisfy many of the properties of the known class of bounded and 
integrable weak solutions. 
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Section |4] studies the actual positivity of nonnegative solutions via quantitative lower estimates for the 



good fast diffusion equation. Precise local lower bounds are contained in Theorem 4.1 The behaviour 



as |x| — )■ oo (so-called tail behaviour) is studied in Section 4.1, and global spatial lower bounds are 
derived as a consequence in Section |4.2[ The merit of the estimates is that they are quantitative and 
most of the exponents are sharp. The lower estimates of this section can be adapted for the exponent 
m = rric separating both fast diffusion subranges, but only when uq S L^^^ for some p > 1. However, we 
refrain from doing this particular case in the present paper since the proof uses some other techniques 
that would lengthen the text. 

The very fast diffusion range < m < ttIc is studied in Section [sj The weighted estimates of 



Theorem 2.2 continue to hold, but this does not allow to obtain the same type of quantitative lower 
bounds since the technique used in the good fast diffusion range does not work anymore. There are 
two problems: on the one hand the L^-L°° smoothing effect does not hold for general initial data, 
on the other hand the presence of the extinction phenomenon makes things more complicated, and the 



extinction time enters directly the estimates of Theorem 5.1 These difficulties have already appeared 
in the standard FDE, s = 1, and were treated in our paper [6]. However the technique used in that 
paper does not extend to < s < 1 and we present here a technique that is based on the careful use 
of weight factors, and in the limit s = 1 gives a simpler proof of the result of |6]. We also study the 
problem of characterizing the finite time extinction in terms of the initial data; thus, we determine a 



class of initial data that produces solutions that extinguish in finite time, see Proposition 5.3, as well 



as a roughly complementary class of initial data for which the solution exists and is positive globally 



in time, see Corollary 5.2 



Section [6] is devoted to study similar questions for the porous medium case. Theorem 6.1 establishes 



local lower bounds of the Aronson-Caffarelli type for all < s < 1. The question of optimal decay as 
l^l — 7- cxD is an open problem; for selfsimilar solutions it is solved in [18J. 

In Section [7] we address a different question that complements our previous results, i.e., the question 
of existence and uniqueness of an initial trace for nonnegative weak solutions defined in a strip Qt = 



(0,T] X Mr. The main results are stated in Theorems 7.2 and 7.3 This result can be combined in 
the reverse direction with the existence of solutions with initial data a nonnegative Radon measure. 
Theorem 4.1 of j l8j . In the final appendix we collect the definitions of weak, very weak and strong 
solutions, together with a number of technical results. 

We still need to mention the relation of these estimates with the linear fractional heat equation, case 
m = 1, for the sake of completeness. The lower bound of Section [6] for m > 1 passes to the limit m ^ 1 
and solves the first case, and this coincides with the limit m f 1 of a part of the estimate for m < 1 
obtained in Section |4] for m < 1. See Proposition |4.3[ 

Notations. Throughout the paper, we fix rric = {d — 2s)/d, mi = d/{d + 2s), pc = d{l — m)/2, 
and "d := l/[2s — d{l — m)] , which is positive if m > rric ■ We will call s-Laplacian of / the function 
— (—A)**/. This is consistent with the use in the standard case s = 1. 



2 Weighted L estimates in the fast Diffusion range 

We will derive weighted estimates which also hold for the standard FDE (i.e., the limit case s = 1). 
When s < 1 the equation is nonlocal, therefore we cannot expect purely local estimates to hold. Indeed 
we will obtain estimates in weighted spaces if the weight satisfies certain decay conditions at infinity. 
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We present first a technical lemma which will be used several times in the rest of the paper. 



Lemma 2.1 Let if G C^(M'^) and positive real function that is radially symmetric and decreasing in 
\x\ > 1. Assume also that ip{x) < \x\~°' and that \D'^(f{x)\ < co|x|~"~^ , for some positive constant a 
and for \x\ large enough. Then, for all \x\ > |xo| >> 1 we have 



|(-A)V(x)| < 



, if a < d . 



C2 log \x 

\x^ 

_C3_ 



dW' 'f(^ = d, (2.1) 
if a > d , 



with positive constants ci,C2,C3 > that depend only on a,s,d and W^Wc'^iM.''-)- ct > d the reverse 
estimate holds from below if (p > 0: \{—Ayip{x)\ > C4|x|~^'^"'"^^'' for all \x\ > \xq\ >> 1 . 



The proof is easy but technical, and is given in Appendix |8.4| for the reader's convenience. We point 
out that the large-decay case a > d is what makes the estimate in the fractional Laplacian case very 
different from the usual Laplacian case. In particular, the s-Laplacian of a nonnegative smooth function 
with compact support is strictly positive outside of the support and has a certain decay at infinity, 
indeed the minimal decay obtained for the {—Ay if when > is compactly supported, 

cf. [12]. A suitable particular choice is the function ip defined for a > as (p{x) = 1 for |x| < 1 and 

f{x) = ^ ^, if|a;|>l. (2.2) 

(l + (|x|2-l)T/8 

We are now ready to present the weighted estimates. 



Theorem 2.2 (Weighted L^ estimates) Let u > v be two ordered solutions to the equation (1.1), 
with < m < 1. Let pr{x) = ip{x/R) where R > and ip is as in the previous lemma with < (p{x) < 
|x|~° for \x\ » 1 and 

2s , 2s 

d < a < d-\ . 

1 — m m 

Then, for all < s,t < oo we have 

)l—m / /. \ 1— m (J If s\ 

< ij^^{u{s,x) -v{s,x))ipR{x)dx\ + j^2l~d{l^m) ^^'^^ 

with Ci > that depends only on a,m,d . 

It is remarkable that the estimate holds for (very) weak solutions, maybe changing sign. Also, it is 
worth pointing out that the estimate holds both for s < t and for s > t. In the limit s — )> 1 we recover 
the well known L^ local estimates for the standard FDE. 
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Proof. • Step 1. A differential inequality for the weighted L^-norm. If -0 is a smooth and sufficiently 
decaying function we have 



_d_ 
dt 



{u{t, x) — v{t, x))ip{x) dx 



--(a) 



<W2 



<(c) 2 



ii-AYu"' -{-Ayv"')^dx 
(n'"-t;'")(-A)^Vdx 
{u-vr\{-Ayij\ dx 

|(-A)^V| 



l—m 



{u — v)il) dx 



i_ 

l—m 



l—m 



dx 



Notice that in (a) we have used the fact that (—A)* is a symmetric operator, while in (6) we have 
used that (ti™ — v"^) < 2^~"^{u — f )™, where = |n|™~^n as mentioned. In (c) we have used Holder 
inequality with conjugate exponents 1/m > 1 and 1/(1 — m). If the last integral factor is bounded, 
then we get 



d 
~dt 



[u{t, x) — v{t, x))tl){x) dx 



[u{t, x) — v{t, x))tlj{x) dx 



Integrating the above differential inequality on (s, t) with s,t > we obtain: 



/ {u{t,x) - v{t,x))ij{ 



x) dx 



l—m 



l—m 



(u(s, x) - v{s, x)) V'(x) dx j < (1 - "l)CA ™ |t - S 



which is (2.3) once we estimate the constant C^, for a convenient choice of test function. 

Step 2. Estimating the constant C^. Choose tjj{x) = (Pr{x) := ip{x/R) = (p{y) , with ip as in Lemma 
"and y = x/R, so that {-Ayip{x) = (-A)V/?(a;) = i?-^'(-A) 



2.1 



K 1-™ 



|(-A)V(y)|i^ 
B2 (^(y)!-™ 



dy + 



|(-A)V(j/)|T^ 

m 

M'* ip{y) 1-'" 
|(-A)V(y)|i^ 

m 

BJr ip{y)^-"^ 



dy 



dy 



where it is easy to check that the first integral is bounded, since (p > k2 > on i?2 , and when \y\ > |xo| 

k3 



with |xo| >> 1 we know by estimates (2.1) that 



-A)V(2/)I 



1 

l — m 



< < 



I l"+T^ ' 
\y\ l — m 

h log 1 2/1 



\y 



if a < d , 



if a = d , 



(2.4) 



d-\-2s — am 

\y\ 



, if a > d , 



therefore ki is finite whenever d — < a < d + ^ . Note that all the constants ki depend only on 
a,m,d. □ 
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Remark. The estimate implies the conservation of mass when (d — 2s) /d = rric < m < 1, by letting 
R — )• oo. On the other hand, when < m < rric solutions corresponding to uq E L^(M'^) n LP(M'^) with 
P > d{l — m)/2s , extinguish in finite time T > , (see e.g. pJJ); the above estimates provide a lower 
bound for the extinction time in such a case, just by letting s = T and t = in the above estimates: 

Moreover, if the initial datum uq is such that the limit as i? — )• +oo of the right-hand side diverges to 
+00, then the corresponding solution n(t, x) exists (and is positive) globally in time, as explained in 



Corollary 5.2 



3 Existence of solutions in weighted L^-spaces 



2.2 



with decay at 



Theorem 3.1 Let < m < 1 and let uq G L^(M"', dx), where (p is as in Theorem 
infinity d — [2s/(l — m)] < a < d + {2s/m). Then there exists a very weak solution u(t, •) E 

(M.'^ , if dx) to equation (1.1) on [0,T] x W^, in the sense that 
rT 



Jm.'^ 



u{t, x)ipt{t, x) dx dt 



Jr"^ 



u"'{t,x){-Ayij{t,x)dxdt, for all ^ e C^{[0,T] x R"^) . 



This solution is continuous in the weighted space, u G C([0,T] : L^(R'^, y^dx)) . 



Proof. Let (p = (pR be as in Theorem 2.2 with the decay at infinity Let < -uo.n £ L^(]R'^) n 

L°°(M'^) be a non-decreasing sequence of initial data uo,n-i ^ uo,n: converging monotonically to uq G 
L-'^(M'^, 9? dx) , i. e., such that j^d{uQ — Un,o)v dx — )• as n — )• oo. Consider the unique solutions n„(t, x) 
of equation (1.1) with initial data uo,n- By the comparison results of |12] we know that they form a 
monotone sequence. The weighted estimates (2.3) show that the sequence is bounded in (R'^ , ip dx) 
uniformly in t G [0,7"]. By the monotone convergence theorem in L^iW^^ip dx), we know that the 
solutions Un{t,x) converge monotonically as n — )• oo to a function u{t,x) G L°°((0,T) : [M.'^ , ip dx)) . 
Indeed, the weighted estimates (2.3) show that when uq G {M.'^ , ip dx) then 

1— m 



u{t, x)p{x) dx 



1—m 



lim 



< lim 



1—m 



Un{t, x)ip{x) dx 

(nn(0,x))(p(x)dxj + CiR'^^^-'^^-^H 
'Uo(x)v?(x)dx ) + CiR'^^^-'^^-'^H 



(3.1) 



At this point we need to show that the function u(t, x) constructed as above is a very weak solution to 
equation ( 1.1 ) on [0, T] x M'^ , more precisely we have to show that for all ^ G ([0, T] x W^) we have 



u{t, x)tpt{t, x) dx dt 



u"'(t,x)(-A)"V(t,x)dxdt. 



(3.2) 



/O JR'i 

By the results of [12J we know that each u„ is a bounded strong solutions, since the initial data 
uo G Li(M'^) n L~(M'^) , therefore for all 4) G C^([0,T] x R'^) we have 



Un{t, x)7pt{t, x) dx dt 



T 



u'^{t,x){-A)'^{t,x)dx dt. 



(3.3) 
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Now, for any ip E C^([0, T] x R'^) we easily have that 



hm 

n— >oo 



J]R<* 



Un{t,x)llJt{t,x)dx 



JR** 



u{t, x)il)t{t, x) dx 



since "0 is compactly supported and we already know that Un(t,x) — t- u{t,x) in i^ioj.- On the other 
hand, for any -0 £ C'~([0, T] x R"') we have that 



lim 

since Un < u and 



T 



<(t,a;)(-A)Xt,x)da;di 



JR'* 



n"(t,x)(-A)^V'(i,a^)da; dt 



< 



JM'' 
T 



< 



< 



x) - x))(-A)"V(i, 2;) dx dt 
|n(t, x) - x)r V7"(x)^^^^3^$^ dx dt 



(^"^(x) 



T 



XJR'i 
T 



|u(t, x) — Un{t^ x)| (/7(x) dx 



|(-A)Xt,x)| 



(p{x) 



1_ 

1 — m 



1— m, 



dx dt 



< c 



{u{t, x) — Un{t, x))ip dx dt — 

/o Jk.'^ 

where we have used Holder inequality with conjugate exponents 1/m and 1/(1 — m), and we notice 
that 

|(-A)^V'(i,^) 



99(x) 



1 

1 — m 



1— m 



dx dt 



< C 



since ■0 is compactly supported, therefore by Lemma 2.1 we know that |(— A)*'(/;(t, x)| < C3|x| ('^+2*)^ 
and the quotient 



-A)^V(i,^) 



1 

1 — m 



< 



C3 

X 1 — m 



is integrable when > d that is when a < d + {2s /m). In the last step we already know that 

Jjgrf(ii(t, x) —Un{t, x))ip dx — )• when ip is as above, i.e. as in Theorem 
in (3.3) and obtain (3.2) . 



2.2 



Therefore we can let n — )■ 00 



For the solutions constructed above, the weighted estimates ( 2.3 ) show that when < G (M*^, dx) 
imply 



/ 



ii(t, x)(y9ij,(x) dx — / u{s,x)ipji{x) dx 



< 2^CiR'^~^^ |t - s\~- 



(3.4) 



which gives the continuity in L^(M'^, c/jdx) . Therefore, the initial trace of this solution is given by 
uq £V^{W^,ipdx) . □ 

Remark. The solutions constructed above only need to be integrable with respect to the weight 
99, which has a tail of order less than d + 2s /m. Therefore, we have proved existence of solutions 
corresponding to initial data uq that can grow at infinity as for any e > . Note that for 

the linear case m = 1 this exponent is optimal in view of the representation of solutions in terms of 
the fundamental solution, but this does not seem to be the case for m < 1. 
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Theorem 3.2 (Uniqueness) The solution constructed in Theorem 3.1 by approximation from below 
is unique. We call it the minimal solution. In this class of solutions the standard comparison result 
holds, and also the estimates of Theorem 2.2 . 



Proof. We keep the notations of the proof of Theorem 3.1 Assume that there exist another se- 
quence < fo,fc G L^(M'^) which is monotonicahy non-decreasing and converges monotonically to 
uq £ L\R'^,^pdx). By the same considerations as in the proof of Theorem 
exists a solution v{t, x) G C([0, T] : L^{W^, tpdx)). We want to show that u ■ 



3.1 



we can show that there 
V, where u is the solution 

constructed in the same way from the sequence no,n- We will prove equality by proving that v < u and 
then that u < v. To prove that v < u we use the estimates 



[Vk{t,x)-Un{t,x)] dx< / [Vk{0,x) - Un{0,x)] 



dx 



(3.5) 



which hold for any Un{t, •),Vk{t, ■ 
get that 



G L"'^(M'^), see Theorem 6.2 of [T2] for a proof. Letting n — oo we 



lim 



[vk{t,x)-Un{t,x)] dx< lim / [^^(0, x) - u„(0, x)] dx = / [vk{0,x) - uo{x)] 

n-5-oo Jjgd -1- J^d 



dx = 



since Vk{0,x) < uq by construction. Therefore also Vk{t,x) < u{t,x) for i > 0, so that in the limit 
A; — cx) we obtain v{t, x) < u(t, x) . The inequality u < v can be obtained simply by switching the roles 
of Ufi and . The validity of estimates of Theorem 2.2 is guaranteed by the above limiting process. 
The comparison holds by taking the limits in inequality (3.5), as it has been done for L^-solutions in 

m- □ 



4 Good fast diffusion range 

The first result of the section will be the existence of local lower bounds. In the proof we will use 
Lemma |8.6[ which is a simple optimization lemma that we state in Appendix |8.5| . We recall that 
rUc '■= d/{d — 2s) and "d := l/[2s — d{\ — m)] which is positive for m > rUc . 



Theorem 4.1 (Local lovi^er bounds) Let Rq > 0, m-c < m < 1 and let < uq G {M.'^ , ip dx) , 
where (p is as in Theorem 2.2 with decay at infinity |x|"°, d — [2s/(l — m)] < a < d + {2s/m). Let 



u{t,-) G h^{W^,ipdx) be a very weak solution to Equation (ll.ip corresponding to the initial datum uq. 
Then there exists a time 

(4.1) 



4. r< R2s-d(l-m) II ||l-m 



such that 



and 



2s 

inf u{t, x) > Ki i?o t^ if < t < t^ , 



flo/2 



inf u{t, x) > K2 



if t > . 



(4.2) 



(4.3) 



The positive constants C*, iCi, K2 depend only on m, s and d> 1. 
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Figure 1: Black: Lower bounds in the two time ranges. Blue: Upper bounds (smoothing effects), which 
has the same behaviour when t > t^ . 



Remarks, (i) The lower estimate for small times is an absolute bound in the sense that it does not 

depend on the initial data (though does depend). 

(ii) We obtain the following expressions for Ki and K2 and C*: 

K2 



K2 :-- 



2s 



d{l — m) 



1 — 71 



1 

1 — m 



and 



d{l-m) {2s^f^^-'^^^-l 



2s 



(2si?) 



d{l-m)-a 



2si? 
1 — m 



a — d 



(4.4) 



2{a-d) + l ojd^'^Cf 



C, = 2s^(u;d2'^I^ 



where Ci > is the constant in the L^-weighted estimates of Proposition 2^ that depends on a, m, d, 
with d < a < d+^, and 1^0 > is the constant in the smoothing effects (4.6), cf. Theorem 2.2 of jl2] . 
(iii) We can always choose a = d/m<d + 2s/m, since 2s > d{l — m). 

Proof. The proof is divided in several steps. 

• Step 1. Reduction. By the comparison principle that it is sufficient to prove lower bounds for 
solutions u to the following reduced problem: 

dtu + i-AYiu"'') = 0, in(0,oo)x™^ 
n(0, •) = uqXBr,^ =uo 



m 



(4.5) 



where rUc <m<l,0<s<l, and iio > . We only assume that < uq € L^(i?/jg) , which implies 
that uq G L^(M°') since supp(uo) ^ Br^ and also that ||iio|lLi(R<*) = ll'"o||Li(Bfl(j) ■ It is i^ot restrictive to 
assume that the ball B^^ is centered at the origin. 
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• Step 2. Smoothing effects. In [12j there are the global — L°° smoothing effects which provide 
global upper bounds for solutions to the Cauchy problem |1.1| . We apply such smoothing effects to 
solutions to our reduced Problem 4.5 to get 



I|2si? _ -"oo II ||2si9 
iLi(lRd) - ^d,?II^OllLi(SflQ) 



II^(*)IIl°°(R'*) < ^IK 

where ■& = l/[2s — d{l — m)] and the constant /qo only depends on d, s,m . 

• Step 3. Aleksandrov principle. We recall Theorem 11.2 of [18], we have that 

u{t, 0) > u{t, x) , for ah t > and |x| > 2Ro . 

Therefore one has that 

h(*)llL-(Md\i32fl„) = _SUp U{t,x) < U{t,0) . 



(4.6) 



(4.7) 



• Step 4. Lower estimates for the L^-norm on an annulus. We combine the L-'^-weighted estimates 
of Theorem 2.2 with the smoothing effects of Step 2: estimates (2.3) read in this context 

/ \ 1— m 



< 



RO 



1— m 



uofR{x)dx] < / u{t,x)^R{x)dx] + CiR'^^^-'^^-'^H (4. 



1— m 



we have chosen R > 2Rq > and ipr{x) = (p{x/R) with (p as in Lemma 2.1 (with the explicit 
form given in formula (2.2)) , so that '■Pr{x) = 1 on Br and < 'Pr{x) < for >> R with 

d — 2s / {1 — m) < a < d + 2s/m , and we recall that Ci > depends only on a,m,d. 



l^o|lL7(^^^)-Cii2'^(i-"^)-2^t< ( / n{t,x)ipR{x)dx 



l—m 



< 



)l—m / 
+ / 



l—m 



(4.9) 



u{t, x)(Pr{x) dx 



(I) + (11) 



We first estimate (I), to this end we observe that if we choose d < a < d + 2s / m we have that 



/ (PRix) dx = fnix) dx + 



Br\B2Rq 



ipR^x) dx 



(fRix) dx + 



1 dx 



BR\B2Rg 



< 



R^^Bn [l + (|x/i?|2-l)4]"/ 



dx + UdR 



+ 00 



1 [1 + (r2 - 1)4" 



a/8 



dr + LOdR" 



<(a) ^dR 



1 [l + (r2_l)4]°/8 

/+00 



dr 



^d-l 



[1 + (r2 - 1)4] 



a/8 



r'^-i-^dr 



OOdR" 



1 + 4'* + 



dr + 1 



4^/8 1 
a-d 4"- 



< .di^'^^-'^-^'R^ 
a — d 



(4.10) 
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where we have used that R > 2Rq and in (a) we have used the fact that ipn < 1 and that 1 + (r^ — l)^ > 
, if r > 4 . Therefore we have 



(I) 



1 

1 — m 



u{t,x)ipR{x)dx <\\u{t)\\i^oor^d\B2R,,) / ^R{x)dx 



< Wd4 



rf2(a-d) + l . 



a — d 



,2ia-d) + l . 



a — d 



where in the last step we have used inequahty (4.7) of Step 3, derived from Aleksandrov principle. 
We now estimate (11) as follows: 



(11) 



1 

1 — m 



u{t, x)(pr{x) dx < \\u 



J Bo n,, 



x) dx 



<(a) ^d^"^ Ro\\u{t)\\i^ao(^d) <({,) UJd^'^RQ^ IMllfiBji^) 



where in (a) we have used that ^r{x) = 1 on Br, 2Rq < R and \Br\ = UdR'^ ■ In (b) we have used the 



smoothing effect (4.6) . Plugging the above estimates into (4.9) gives 

id Tjd i 



I , 111— m /-( r>d(l—m)—2s 4. ^ 



^d 2^ RqIoo II'"o|Ili'(b^^ 



1— m 



+ 



^d^ 



^d(l-m)i? 

d'2{ a-d) + l 
a — d 



l—m 



(4.11) 



or equivalently 



l—m ■ 



1 



(4.12) 



< 



Wd4' 



,2(a - d) + 1 
a — d 



l—m 



ni-™(t,0) 



Step 5. Optimization. The previous estimate (4.12) is useful only if we can make sure that the 



left-hand side has a positive lower bound. Let us write inequality (4.12) as 

A{t) Bt 



with 
A(t) 



M 



F{t,R) :-- 



C 

^d(l-m)i? 



' Rdil-m) r2s 



< 



,2{a -d) + l 
a — d 



l—m 



u^-"(t,0), 



LOd 2'^RqIoo \\uo\\l4br^^) 



l—m 



(4.13) 



B = Ci (4.14) 



where Ci > is the constant of L^-weighted estimates of Theorem 2.2, and I^o > is the constant of 
the smoothing effects (4.6) of Step 2. We now optimize the function F as in Lemma 8.6 so that there 
exists 



2519 



(J \ d{l-m)7j 

M 



2s^(ud2^I^]"" R^ llnoll^T 



l—m 



(4.15) 
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and 



2sBt 



d{l-m)A{t) 

2s (2si?)2^'^ 



2s 



d{l - m) (2si?)'^(i-™) - 1 



'^^olooCf >(a) 2Rq 



(4.16) 



d{l - m) (2si?)'^(i-'") - 1 

where in (a) we have used that the constants loo > and Ci > are constants in the upper bounds 
(4.6) and (4.8) respectively, so that we can chose them to be arbitrarily large to fulfill the condition 
R{t*) > 2R. Therefore for all t > t=K we have that 



(2(a -d) + l 
a — d 



1— m 



u'-"'{t,0) > F{R{t),t) 



2s 



d{\ — m) 



d{l — m) 



2s 



2si? 



Ait) 



2si? 



> 



2s 



1 



^d{l — m) 

since A{t) > A{t.t) for all t>t^, and it is easy to check that 



A{U) 



1 



{2s'd) 



d(l-m)'& 



1 1—m 



d{\ — m) 
2s 



A{U 



1 



^d(l-m)i? ^d(l-m)i? 



since we recall that 2s'& > 1. Summing up we have obtained 



u{t, 0) > K2 



U, l|2si? 



(4.17) 



for all i > f,, > , where K2 only depends from a, m, s, d and takes the form 



K2 :-- 



2s 



d{l — m) 



1 

1 — m 



d{l-m) (2si?)^(^-'")^^-l 



2s 



{2sd) 



d{l-m)'d 



2s-d 
1 — m 



a — d 



2(a - d) + 1 ujd^dcf 



(4.18) 



Note that in the limit m — >• 1 the constant K2 — )• 0. By a standard argument it is easy to pass from 
the center to the infimum on -6^^/2(0) ™ the above estimates. 

• Step 6. Positivity backward in time. Using Benilan-Crandall estimates which depend only by the 
homogeneity of the equations, cf. ^ 



ut < 



u 



1 — m)t 



(4.19) 



we can prove positivity in the time interval [0,t=K]. These estimates in the fractional case has been 
proven in [12], and imply that the function: u{t, x)t~^^^^~'^^ is non-increasing in time, thus for any 
t G (0, t^:) and x G ^_Ro/2(0) > inequality ( 4.17[ ) gives 



2sS 



i{t,x) > 1 ti-" > K2 — ^ 1 ti-" 

J. 1 — m J. 1 — m 



K2 



2^+^sd iuJdloo)^ 



den 
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where K2 > is given in (4.18) , and is given by (4.15), and it is easy to check that 

1 



U, ||2si9 
F0|lLi(i3„ ) 



u, \\2s-a 

F0|Il1(B«J 



The proof is concluded. □ 



Remark. This lower estimate holds in the limit m — )■ 1 and gives lower estimates for the linear 
fractional heat equation of the form 



Proposition 4.2 Let u > be a weak solution to the Cauchy Problem (1.1), corresponding to uq G 
L-'^(M'^) and m = 1. Then i9 = l/2s > and the estimate says that for given Rq and t^: := Ci Rq^ , then 



inf u{t, x) > K2 



Ro/2 



l^ollL^Bflp) 
td/2s 



if t > . 



(4.20) 



The positive constant K2 depends only on Ci, s and d. 



The proof is easily obtained from the integral representation of the solution. 



4.1 Minimal space-like tail behaviour 



As a corollary of the previous lower bound, we obtain a quantitative bound from below for the space- 
like behaviour of any nonnegative solution. We consider a solution that has a certain initial mass M in 



the ball of radius 1 and apply the result of Theorem 4.1 after displacing the origin of space coordinates 
to a point xq with |xo| >> 1- We then consider the formula (4.1) for the critical time with center xq 
and radius Rq = |xo| +2, so that the ball Bii^^{xq) contains the mass M mentioned above. As Rq 
also t* — )• 00. We can therefore use the lower bound (4.2) to get an estimate of the form 



00 



u(t,xo) > G{uQ,t)\xQ\-'^'/^^-'^^ 



(4.21) 



where G{uQ,t) is given in (4.2). According to the results of |I8] the Barenblatt solutions have this 
precise spatial behaviour in the range rric < m < mi, with mi = d/{d + 2s), therefore the asymptotic 
estimate is sharp in this range. 



4.2 Global spatial lower bounds in the case mi < m < 1 

We would like to prove that the solution can always be bounded from below by a Barenblatt solution, 

gives 
In 



so the lower bound will be sharp. In the range ruc < m < mi the lower bound of Theorem 4.1 



sharp lower bounds with the same tails as the Barenblatt solutions, as explained in Section 4.1 



the range mi < m < 1 the lower bound given by (|4.21|) is not sharp and the following Theorem 4.3 
(respectively Proposition 



4.2 



when m = 1) proves that any solution with data in iJ'iJ 



can always be 

bounded from below by a Barenblatt solution (respectively by the fundamental solution when m = 1). 
See Figure [2] 
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Figure 2: Lower bounds for the spatial decay rates of solutions. Recall that rric = {d — 2s) /d and 
mi = d/{d + 2s) 



Theorem 4.3 (Global Lower Bounds when mi < m < 1) Under the conditions of Theorem 4.1 
we have in the range nii < m < 1 

u{t, x) > ^ when\x\»l. (4.22) 

valid for all < t < T with some bounded function C > that depends on t, T and on the data. 

Proof. The proof consists of several steps. 

• Step 1. We begin under the extra assumption that uq{x) > 2c > in a ball, that can be taken to 
be -Bi(O) by scaling. Therefore, there exists a ti such that u{t,x) > c for all < t < ti and all |x| < 1. 
We also assume that uq is continuous and goes to zero uniformly as |x| — )• oo. 

Consider the function uo^s{x) = uq{x) + e, and let Ui.{t,x) > u{t,x) be the corresponding solutions. 
By the usual theory, [12j, we know that Ue > e, Ue — e ^ L^(]R'^), since no G L^(M'^). Moreover, it is 
proved in the theory that — )• e as — )• oo for every t > 0. 

• Step 2. Let u = u*{t + T,x), where u* is the Barenblatt solution with mass M > 0. We refer to [18] 
for a complete discussion about Barenblatt solutions. By choosing the mass M := udx very small, 
we can find r = r(e) > so that u{t, x) < c/4 for \x\ > 6 and < t < ti, and u < e/2 for \x\ > 1. 

• Step 3. We compare both continuous solutions in the exterior domain Q = {x : \x\ > 1}. At the first 
time where u touches Ue from below at a point |x| > 1, we have dt{us — u) < 0. Let now w = — u"^ 
to get 



/ A N« / N , f w(x) — w(y) - 
jRd \x — y\ ^ 



, w(x) — win) , , f w(x) — win) , 

{\x\<i} \x-y\'^+^' J{\x\>i} \x-y\'^+^^ 
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We want to prove now that both /i < and I2 < 0, which leads to a contradiction. In this way we 
conclude that u < Us for all < t < ti and all |x| > 1. 

The fact that I2 < comes from the fact that ■w{x) = by our choice of x and w{y) > since (x,t) 
is the first contact point. Due to the fact that w{y) > near \x\ = 1 and w is continous we get I2 < 0. 

As for /i, the denominator is like a constant in that domain and we have to estimate w{y). We know 
that for 6 < \x\ < 1 we have Ue > c and u < c/4, hence w{y) > C c™" > and this contributes to the 
integral something that is like —Cc"^, which is not small. In the small ball |x| < 5 we use the worst 
case estimate —w{y) < u and has mass at most M which is small, this contributes at most a 

bad term of order 

C [ < CM™^'^^!-™), 

J\x\<S 

which is small if 6 and M are small (here we use m < 1). Therefore Ii < 0. 
Moreover, one has to ensure that u{0,x) < Ue(0,a;) for > 1. Since > e and 

at least for sufficiently small r, recall Step 2. By the parabolic comparison theorem we conclude that 
u<Ue for all < t < ti and all > 1. 

• Step 4. We finally let e — t- and also r may go to zero, and we obtain that u*{x, t) := lime_j.o = u, 
therefore we can conclude that u[x,t) > c/\x\'^~^'^'^ when \x\ » 1 and t = ti. 

• Step 5. Once we have obtained the spatial lower bound at times t < ti, then we can compare with 
a Barenblatt solution and continue the lower bound for all times, to finally get that the spatial tail of 
the solution u can be bounded from below by u > c/|a;|°'+^* when |x| >> 1. □ 



5 Very fast diffusion range 



2.2 



continue to 



In the very fast diffusion range < m < rric, the weighted estimates of Theorem 
hold, but this does not allow to obtain quantitative lower bounds since technique used in the good fast 
diffusion range does not work anymore. One problem is that the smoothing efi^ect does not hold for 
general initial data, therefore the optimization of Lemma 



8.6 



is no more valid, since 2s < d{l 



m) 



in this range. Hence the need for new weighted estimates, in the form given in Step 3 of the proof of 



Theorem 5.1 below. Another problem typical of this range of exponent is the presence of the extinction 



time, which enters directly in the estimates of Theorem 5.1 We present here a technique that is based 
on the careful use of weight factors. 



Theorem 5.1 (Local lower bounds I) Letu be a weak solution to the equation (1.1), corresponding 
to uo G Li(M'^) n LP-{R'^) with < m < = d/{d - 2s), Q < s < I and let pc = d{l - m)/{2s). Let 
also T = T{uo) be the finite extinction time for u. Then for every Rq > 0, there exists a time 



such that 



^ p2s-d(l-m)|| ||l_m 



inf u{t, x) > K- 



Ro/2 



Rn 



2^m(l-m) 



if 0<t<t^, 



(5.1) 



(5.2) 
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where C* and K are explicit positive universal constants, that depend only on m, s, d. 



The expression of the constants is 



C ■= ^"'-^^^^ K ■= ( ^ ] (5.3) 

^* ■ 4d+l-2s ' V 43'^+l-2s d) ' 

where kg^d is the constant of the representation formula (p{x) = ksA J^d \^J'yp-2s and ujd is the volume 
of the unit bah. 



Proof, of Theorem 5.1 It is divided into several steps as follows. 

• Step 1. Reduction. By the comparison principle that it is sufficient to prove lower bounds for 
solutions u to the following reduced problem: 

dtu + {-A)%u'^) = 0, in(0,oo)xM^ 
u{0,-) = uqXBr^ = uo, inM'^, 

where m>l,0<s<l, and Rq > . We only assume that < uq £ L^(i?/j„) , which implies that 
Up G L^(M'^) since supp(uo) Q and also that ||?^||li(R'*) = ll^ollLi(Bflg) • It is not restrictive to 
assume that the ball iJ/jg is centered at the origin. We call Mq = \\uo\\i^i(^Biig)- 

• Step 2. Aleksandrov principle. We recall Theorem 11.2 of jl8j . In view of the fact that the initial 
function is supported in the ball Bjig{0), we have that 

u{t, 0) > u{t, x) , for all t > and |x| > 2Ro . 

Therefore, one has 

sup u{t,x) < u{t,0) . (5.5) 



• Step 3. L-*^ Weighted estimates. Choose a test function > such that —{—A)^ip = p with /? = 
on B2Rq and on and < /) < 1 on the annulus A := Br^ \ i?2i?o' with < 2Rq < Ri, and Rq as in 
Step 1, such that supp(no) ^ -B/j,,. Using the explicit representation of (/? in terms p and the integral 
kernel K{x,y) = ks^d\x — we get the estimates 

^(-) = I ^^r^ > > «=o > , for al. x 6 S«.(0) . 

since |x — y| < Rq + Ri. We can always choose p > 1/2 on the smaller annulus Aq = -B2i?,o+3(-Ri-2_Ro)/4\ 
^2Ro+(Ri-2Ro)/4 ^ ^, so that 

IIpIIi = / P{x) dx> [ p{x) dx > = ^|52fio+3(iJi-2iJo)/4 \ B2Ro+(Ri-2Ro)/a\ 

J Al JAo ^ ^ 



2 



2i?o + J(^i - 2i?o) + iRx - 2Ro)] - (2R0 + \{Ri- 2Ro ' 



>'f{R,-2R^r 



since (a + b) — a > b for any a,b > 0. Then ko > takes the form 



, ks.dOJd {Ri-2Roy 

^''■= 2 {R, + Ror-'^ ^'-'^ 
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Now we observe that letting T = T(uq) > be the finite extinction time for the reduced problem (5.4), 
we obtain 

UQ{x)ipdx= / u{0,x)ipdx— / u(T,x)ipdx = — / / dru{T,x)(p{x) dxdr 

JRd. J^d Jq J^d 

= r [ {-Ay{u"'{T,x))^{x)dxdT= r [ n™(T,x)(-A)V(x)dxdr 

Jo JR'i Jo JRd. 

rT f 

u^{t, x)p{x) dxdr 



J A 

' u'^ir, x)p{x) dxdr + I I n"(T, x)p{x) dxdr := (/) + (//) 



(5.7) 



10 J A Jtt. J A 

where Q <t^, <T will be chosen later. 
Next we estimate (I). We first observe that 



u"'{T,x)p{x)dx< n™(r,x)dx<|B/jji-'" / u(r,x)dx 
Jbr^ \Jbr^ ) 

<\Br,\'-"'([ u{T,x)dx) <\Br,\'-"^([ uo{x)dx 

<\Br,\'-^M^ 

since 0<m<l,0<p<l, and in the last step we have used the fact that ||'u(i)llLi(R'*) — II^o||l1( 
which has been proven in [12j , together with the fact that Mq = 11^*01^1(5^^) = II''^o||li(r<*)- Therefore, 

(/) := r* f u"'{T,x)p{x)dxdT KIBr^I^-'^UM^. (5.8) 
Jo J A 

We now estimate (II) by using the Aleksandrov principle: 

(//) := I" I" n"(T, x)p{x) dxdr < [ [ """(r, x) dxdr =(„) (T - U) [ u'"(ri, x) dx 



t, JA Jt, JA 



< (r-t,)|yl|supn"(Ti,x) < (r-i,)|^|u"^(ri,0) 

where in (a) we have used the mean value theorem for the function U{t) = J^u'^{t, x) dx so that 
there exists a ri G [i*,?"] such that C/(r)dr = (T — t*)C/(ri). In (6) we have used the Aleksandrov 



principle, which gives supM'"(ri,x) < n™(ri,0). Summing up, we have obtained, joining (5.7), (5. 
and ([5^ 

/ uo{x)y^dx<\BR,\^-"'UM^ + {T-U)\A\u'^{n,0) (5.10) 

jRd 

for some ri G In addition, we have J^aUo{x)(pdx > Mq/cq- We finally remark that from 

inequality ( 5.10| ) we get a lower bound for the extinction time, just by letting t,,, = T in formula (5.10): 

koMo< / uo{x)ipdx<\BR,\^-'^TM^, that is T > ko ° (5.11) 
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Step 4. Choosing the critical time t^,. We now choose t=K to be small enough, more precisely 



t. 



2 \Br, 



l—m 



1—m 



< T, 



we note that t^, <T follows by (5.11). With this choice of t^,, inequality (5.10) becomes 
^Mo = koMo-\BR,\'-^t,M^ < (r-t,)|^|n"^(ri,0) <T|^|7x"^(ri,0) 

which is the desired positivity estimate at a time ri G [t-t , T] , namely 

/Co Mo 



2r \Br^ \ B2Ro\ 



(5.12) 



(5.13) 



(5.14) 



• Step 5. Positivity backward in time. Using Benilan-Crandall estimates which depend only by the 
homogeneity of the equations, cf. |3] 



ut < 



u 



(5.15) 



;i - m)t 

we can prove positivity in the time interval [0,ri]. These estimates in the fractional case has been 
proven in [T^, and imply that the function: u{t, x)t~^^^^~"^'^ is non-increasing in time, thus for any 
t £ [0, Ti] we have that 



1 1 

^ 1 — m. ^ 1 ~ m 

u{t,0) > — —u{ti,0) > —u{ti,0) > 



ka 



TT- 



(i?l — 2Rq 



koMo 



2T \Br-^^ \ B2Rq\ 



X 1 

m ^ 1 — n 



(5.16) 



since < ri < T . Moreover we have that 



2^m{l-m) 



n(t,0) > 



k 



s,d 



{Ri - 2Rc 



\d-i 



4(i?i + Ro)d--^' d{2Ro) 



T \ 1/m / o 



Ad 



2Ro 



i-l 

-JL 1 



fl-m 1_ 
2^m(l-m) 



(5.17) 



2^m.(l-m) (i?^ + Rq 



where we have used the numerical inequality a'^ — b'^ < da'^~^{a — b) , valid for any a = Ri > 2Rq = b to 
pass from ( |5.16 ) to (5.17). By a standard argument it is easy to pass from the center to the infimum 
on -B_Ro/2(0) in the above estimates. The proof is concluded once we let Ri = SRq. □ 

Remarks, (i) This result can be written alternatively as saying that there exists a universal constant 
Ki = mayi{K~"^,Cy^^ "^^} such for all solutions in the above class we have: for any < t < T and 
R>0 



R'^ 



< Ki 



+ 



TT-- 



R^- 



inf u™(t,x) 



(5.18) 



This is easy to prove: by the previous Theorem, we have that either < t, that is 



Rd 



1 

l — m 



< 
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or that < t < and (5.2) holds, namely 

ll«o||Li(Bfl) 



< 



inf vJ^{t,x) 



l/(l-m) 



} we get (5.18) . 



therefore, letting Ki = max{i(' C* 

This equivalent version is in complete formal agreement with similar estimates proved by the authors 
in [6], in the case s = 1. However, our proof below differs from the one in [6], and provides an alternative 
proof when s = 1. On the other hand, here we are considering solutions to the Cauchy problem, while 
in [6] we consider local weak solutions (i.e. without specifying boundary conditions). These estimates 
have been called Aronson-Caffarelli estimates in [6j, when s = 1, since they are quite similar to the 
one that can be obtained for m > 1, see Section [6] Finally we shall remark that in Section 5.1 we will 
obtain quantitative upper estimates on the extinction time, and this will help to eliminate T from the 
above lower estimates. 

(ii) By comparison it is easy to prove that this estimates hold for a larger class of solutions, more 



precisely for the class of very weak solutions to the Cauchy Problem (1.1) constructed in Theorem 3.1 
Sectionjsj This implies that the positivity result holds for solutions u {t, •) G L^(M'^, ipdx) corresponding 



to initial data < uq G {R'^ , tp dx) , 
d - [2s/(l -m)]<a<d+ (2s/m). 



where 99 is as in Theorem 2.2 with decay at infinity \x 



Once comparison is used, we can use as T the extinction time of the reduced problem 5.4 in Step 1 
of the above proof. In this way the quantitative result applies to solutions u that may not extinguish 
in finite time. Therefore we can interpret T as the minimal life time for the solution u{t, •), a concept 



that was already introduced by the authors in [6] , for which formula (5.1 ) provides a quantitative lower 
bound, namely 



2s-d{l-m) I 
I 



(5.19) 



2.2 



Corollary 5.2 (Solutions that do not extinguish in finite time) Let < m < nic and consider 
an initial datum < uq £ L-'^(M'^, ipdx), where ip is as in Theorem 
Assume moreover that 

liminf 



+00 . 



in particular, when uq G L^(M'^). 

(5.20) 



PoIIlhBb) 

Then the corresponding solution u{t, x) exists and is positive globally in space and time, hence does not 



extinguish in finite time. Moreover the quantitative lower bounds (5.2) of Theorem 5.1 hold for any 



< t < witht^ given in (5.1) andT = T{uqXBb,q) < +00 is the extinction time of a reduced problem. 



Proof. If we consider an initial data with that behaviour at infinity, then by Theorem |3.1| there exists 
a very weak solution. By letting R 
minimal life time T{uqxBb) 
since < m < rric. □ 



+00 in the above lower bound (5.19) for T, to conclude that the 
00, recalling that in this very fast diffusion range we have 2s < d{l — m), 



Remark. A practical assumption on the initial datum uq that implies (5.20) is 



liminf \x\"^-^uq{x) 

\x\-^+oo 



+00 . 



(5.21) 



In view of Proposition 5.3 below, the exponent is sharp. 
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5.1 Estimating the extinction time. 



We next estimate the extinction time in terms of the initial data, extending a classical result of Benilan 
and Crandall [2]. This is needed to eliminate the dependence on T in the above lower estimates when 
we consider initial data in L-^(]R'^) n U'^iW^). For a detailed study of extinction time in the standard 
fast diffusion equation, see [Kj. 



Proposition 5.3 (Upper bounds for the extinction time) Let u he a weak solution to the equa- 
tion (1.1), corresponding to uq G L-'^(M'^) n LP'=(M'^) with < m < rric = d/{d — 2s), < s < 1 and let 
Pc = d{l — m)/{2s). Then for all < s <t the following estimate holds true 



\u{t, x)\^'' dx 



< 



\u{s, x)\^'' dx 



4m[d{l -m) - 2s\ 



t - s) 



d{d - 2s)Sl 

Moreover, there exists a finite extinction time T > which can be bounded above as follows 

d{d - 2s)S^ 



T < 



4m[d(l -m) - 2s] 



u, 111—™' 



(5.22) 



(5.23) 



Proof. The proof presented below is analogous to the one of Theorem 9.5 of [12] , but here we pay 
attention to the quantitative estimates . We multiply the equation by with p > 1, and integrate 

in M'^. Using Strook-Varopoulos inequality (8.6) in the form (8.7), we get 



d^ 
dt 



{t,x)\Pdx = -p \u\P-^u{-Ay{\ur-^u) dx 



< 



< 



Amp{p — 1) 
[p + m — 1)^ 

4mp{p — 1) 
{p + m — 1)^5^ 

Sobole^ 

make the choice p = Pc = d{l — m)/2s, so that p, 



-A)2 n 2 



dx 



(5.24) 



d-2s 

d(p+m-l) \ d 

\u\ <^-2s dx 



I 1 p+m — 1 

where in the last step we have used the Sobolev inequality (8.8) applied to / = |n| 2 . Now we 



m < nic = d/{d — 2s) , and inequality (5.24) becomes 



^^^d-2s ' know that Pc > 1 if and only if 



d r ^ Am\d{l - m) - 2s\ 

— / ki(i, x) P'= dx < LA- ^—-^ — - 

dt J^i - 2s{d-2s)S^ 



|n(t,x)|P'= dx 



(5.25) 



Integrating the above differential inequality on {s,t) gives both (5.22) and inequality (5.23) . n 

Thanks to the above estimates we can get rid of the extinction time T in the lower estimates of 
Theorem 15.11 



Theorem 5.4 (Local lower bounds II) Let u be a weak solution to the equation (1.1), correspond- 
ing to Uq e L^(M^) n LP'=(M'^) with < m < rric = d/{d - 2s), < s < I and let pc = d{l - m)/{2s). 
Then for every ball ^^fi'"^ exists a time 



^ p2s-d(l-m)|| ||l-m ^n^(,, ^ 



<C\\uo\\lr 



(5.26) 
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where we recall that T{uq) is the finite extinction time, such that 



7;„ rn 1 

"U T 1( fi„ \ f 1-m 

inf u{t,x)>K, 



if 0<t<t^, 



(5.27) 



where C* and Ki are explicit positive universal constants, that depend only on m, s, d. 
The expression of the constants 



Ko ■.= K 



4m[d{l - m) - 2s] 
d{d - 2s)52 



m(l — m) 



C :-- 



d{d - 2s)52 



(5.2S 



4m[d(l -m) - 2s] 

where C* and K are as in (5.2) and kg^d is the constant of the representation formula '^{x) = 
k J f Pfa) rl?j 

Remark. This result can be written alternatively as saying that there exists a universal constant 
= mayi{K2™',Cy^^ "^^} such for all solutions in the above class we have: for any < t < T and 
R>0 



No||li{Bb) 



R'^ [l V ||uo||lpc(b 



1 



inf u"'{t,x) 



(5.29) 



This equivalent version is in complete formal agreement with similar estimates proved by the authors 
in [6], in the case s = 1. 



Proof of Theorem\5.4\ The proof is the same as for Theorem 5.1 , as far as the first 3 steps are concerned. 



At the end of Step 3, we need to bound from above the extinction time for the reduced problem (|5.4|) 

d{d - 2s)5| 



with the estimates (|5.23|) which give 

|iio(x)|^'' dx 



T < 



d{d - 2s)Si 



4m[(i(l -m) - 2s] 



Am[d{l -m) - 2s] 



|iio(x)|^'' dx 



2s 

d 



(5.30) 



since supp(no) C . Then the proof follows simply by replacing T with the above upper bound. □ 



6 The Porous medium case 

Lower estimates for nonnegative solutions of the standard porous medium equation were obtained in 
Aronson-Caffarelli in a famous paper yy. We want to show in this section how such a priori estimates 
extend to the fractional version considered in this paper. 



Theorem 6.1 (Local lower bound) Let u he a weak solution to Equation (1.1), corresponding to 



uq G L^(]R'^). and let m > I. We put "d := l/[2s + d{m — 1)] > 0. Then there exists a time 

I — (m— 1 



such that for every t>t^ we have the lower hound 

inf u{t,x) > K ° ^'^-^^^ (6.2) 
valid for all R > 0. The positive constants C and K depend only on m, s and d, and not on R. 
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Remark. This result can be written alternatively as saying that there exists a universal constant 
Ci = Ci{d, s, m) such for all solutions in the above class we have 



no 



(6.3) 



This equivalent version is in complete formal agreement with Aronson-Caffarelli's estimate for s = 1. 
However, our proof below differs very strongly from the ideas used in Aronson-Caffarelli's case since 
we cannot use the property of finite propagation of solution with compact support, which is false for 
s<l. 

Proof. It is divided into several steps as follows. 

• Step 1. Reduction. By the comparison principle that it is sufficient to prove lower bounds for 
solutions u to the following reduced problem: 



+ (-A)^(u™) = 0, in(0,oo)x 
n(0, •) = uoXBr,, =uo, inW^ , 



(6.4) 



where m>l,0<s<l, and i?o > . We only assume that < € L^(i?/j„) , which implies that 
uq G L^(M'^) since supp(uo) ^ Br^^ and also that ||'Uo||li(R'') = II^o||li(Bh(j) • It is not restrictive to 
assume that the ball S/jy is centered at the origin. 

• Step 2. Smoothing effects. In |12j there are the global L^-L°° smoothing effects, which can be 
applied to solutions to our reduced Problem 6.4 as follows: 



< 



^ll!fO|lLi(IRd) - ^IF0|lLi{iJB„) 



(6.5) 



where = l/[2s + d{m — 1)] and the constant Iqo only depends on d, s,m . 

• Step 3. Aleksandrov principle. We recall Theorem 11.2 of |18j. In view of the fact that the initial 
function is supported in the ball Bfi^{0), we have that 



Therefore, one has 



u{t, 0) > u{t, x) , for ah t > and \x\ > 2Rq . 
sup u{t, x) < u{t, 0) . 



(6.6) 



• Step 4. Weighted estimates. If ■0 is a smooth, nonnegative, and sufficiently decaying function, we 
have 



dt 



u{t, x)ip{x) dx 



((-A)"n™)Vda; 



--(a) 



n"(-A)>dx 



(Rd) / u{t,x)dx 



< 



jrn—l 

oo 



|2si?(m-l) 



-A) 



m— 1 

oo 



A\s I II II ||2si?(m-l)+l 

-A) ^|Il-(R<^) WML^Ba,) 



(Kd) / uo{x)dx 



^di?(m-l) ' 
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Notice that in (a) we have used the fact that (—A)'' is a symmetric operator, while in (b) we have used 



the smoothing effect (6.5) of Step 2 and the the conservation of mass: jj^d x) dx = f^aUo{x)dx, 
for alH > 0, together with the fact that supp(iio) ^ -B/Jq. We refer to | |12| for a proof of the smoothing 
effect and of the conservation of mass. Summing up, 



d 
~dt 



u{t, x)ip{x) dx 



< 



since d'&{m — 1) = 1 — 2s'Q. Integrating the above differential inequality on (0, t) with t > we obtain: 



K[uo,iIj] 
2sd 



t^^^ < I u{t,x)tl'{x)dx- I u{0,x)tl'{x)dx < 



2sd 



We will use this in the form 



u(0, x)^{x) dx - < I u{t,x)i^{x)dx. 



Moreover, if ^ G and Ri > 2Ro, we have 



u{t, x)tP{x) dx 



u{t,x)ip{x) dx + / u{t,x)il}{x) dx 



<(a) \Bji^\ sup u{t,x)+ sup u{t,x) / ■ip{x)dx 
\x\<Ri xm''\B2Ro Jb-j^ 



<{b) I^Ril^ll^^ollLitsflo) / ^ i^{x)dx 



(6.7) 



(6.8) 



where in (a) we have used the fact that < I, Ri > 2Rq , and that G L (M ). In (6) we have 



inequalities (6.7) and (6.8) , we obtain 



used the smoothing effect (6.5) of step 2 and the Aleksandrov principle of Step 3. Putting together 



n(0,.)^(x)dx-^|^t- 



^RiIJ^II^o|Ili(b«^) <^(i,0) / ^ ^{x)dx 



(6.9) 



Next, in order to estimate K[uo,^l^] i^i & convenient way we take ^{x) = (f){\x\/R) with <j) as in Lemma 
, we have |(— A)'^?/;| < cs^Z"^* , for some constant C3 = C3{d, s) . Then, 



2.1 



jm—l 

'-00 



KIuq^iIj] = I[ 
When R> Rq and Ri > 2Rq, we arrive at 



A / II II II— v- . - ^ " c 

A) VIlL-CRd) II^0|Il1(b^^) < 



m— 1 

00 II ||2si3(m-l)+l 
II^0|Il1(b^^) 



(6.10) 



Po||li(b«,) 



II ||2si3{m-l)+l 



2s?? 



t 



2s-a 



-u;,Rf^\\uo\\-f\%^^)<u{t,0) I ^ ^l^{x)dx 
< u{t, 0)R'^ / ip{x) dx = CiR'^u{t, 0) . 



(6.11) 
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• Step 5. Choosing the parameters. We want to choose t > 0, R > Rq and R > 2Rq so that the 
left-hand side of (6.11 ) is larger than ||iio||Li(_B^^)/2, which will then give the desired bound from below 
for u(0,t). We first make the choice 

II ||(i(m-l)i9 

which will satisfy the condition Ri > 2Rq if and only if t > where 

tt = cMuo\\l^~^^J. C5 = 2i+(2/<^(u;,/oo)^/'^. (6.13) 
Now we can make the second choice, R has to be large enough, for instance: 

II ||i?(m-l / 4C3ioo \ (a^A\ 



Both choices will give for t > t^^ the lower bound 

II^oIIlMBho) 



<u(t,0). 



2c4R'^ 

which can be rewritten as 

....... --2c4c: 



C7 ,„ ^" < n(t, 0) , for any t > , with cy = -i . (6.15) 



By a standard argument it is easy to pass from the center to the infimum on -6/20/2(0) the above 
estimates. □ 



Remark. In the limit m — )• 1 of the estimate of Theorem 6.1 we obtain the result of Proposition 4.2 
for m = 1. 

Open Problem. To calculate the positivity of the solutions for small times is not known yet. 



7 Existence and uniqueness of initial traces 



The existence of solutions of the Cauchy Problem (1.1)- (1.2) can be extended to the case where the 
initial datum is a finite and nonnegative Radon measure. We denote by A^+(M'^) the space of such 
measures on W^. Here is the result proved in Theorem 4.1 of |18| . 

Theorem. For every ji G A^"'"(M'^) there exists a nonnegative and continuous weak solution of Equation 
(1.1) in Q = (0, 00) X taking initial data ji in the sense that for every cp G C^(M'^) we have 

n(t,x)^(x)dx = /^(x)dMx). (7.1) 



lim 

t-s>0+ 



In this section we address the reverse problem, i.e., given a solution to find the initial trace. In the 
case s = 1 such question was solved thanks to the works of Aronson-Caffarelli [Ij, Dahlberg-Kenig [S], 
Pierre [13J and others, see a presentation in [T7|, Chapter 13. 
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Lemma 7.1 (Conditions for existence and uniqueness of initial traces) Let m > and let u 

be a solution to equation (1.1) in (0,T] xM"^. Assume that there exist a time < Ti < T, some positive 
constants Ki, K2, a > and a continuous function u : [0, +00) — )• [0, +00), with a;(0) = such that 



sup / u{t, x) dx < Ki 
ie(0,Ti] JBnixo) 



y R> 0, xo G 



as well as 



u{t, x)(p{x) dx 



< 



u{t' , x)ip{x) dx 



+ K2Uj(\t-t'\ 



(7.2) 



(7.3) 



for all < t,t' < Ti and for all ip € C^(M'^) . Then there exists a unique nonnegative Radon measure 
fi as initial trace, that is 



(/9d/x= lim / u{t,x)(p{x)dx , for all (p e {R'^) . 



Moreover the initial trace /i satisfies the bound (7.2) with the same constant, namely fi{Bji{xQ)) < Ki . 

Notice that the constants Ki and K2 may depend on u and ip, usuaUy through some norm. 

Proof. The proof is divided in two steps in which we prove existence and uniqueness of the initial trace 
respectively. 

• Step 1. Existence of the initial trace. Hypothesis (i) easily implies that 



lim sup 



u{t, x) dx < Ki 



V i? > , xo G 



t^0+ JBr{xo) 

Moreover, it implies weak compactness for measures (to be more precise, weak* compactness, see 
Theorem 8.8 in the Appendix 8.6), so that there exists a sequence tfc — )• 0"*" as A: — )• 00 with Q < t^ < Ti , 
and a nonnegative Radon measure /i so that 



lim 



u{tk,x)(p{x) dx 



/ (pdpL for all ip G C°( 



The bound (7.9) on the initial trace: ii{Br{xq)) < Ki follows from the above bound on the lim sup. 

• Step 2. Uniqueness of the initial trace. The initial trace whose existence we have just proved may, 
of course, depend on the sequence t^. We will now show that this is not the case, thanks to hypothesis 
{ii). Assume that there exist two sequences tk — >■ 0"^ and — t- 0^ as A; — >• cxd , so that u{tk) — >• and 
u{t'^) —J- z^,with n,u e M~^{R'^). We will prove that 



(/? dfi 



ipdv for all £ 



(7.4) 



so that fi = f as positive linear functionals on C^(M'^). Then by the Riesz representation theorem (cf. 
Theorem 8.9) we know that fi = u also as Radon measures on R'^. Therefore, it only remains to prove 
( pni ) : hypothesis (ii) implies that for any t,t' >0, with < t + t' <Ti <T, and any (p G C^(M'^) we 
have uj{\{t + t') - t\) = uj{t') and 

u(t,x)ip{x)dx < I u{t + t',x)ip{x)dx + K2Uj{t') . (7.5) 
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First we let t = and t' > to be chosen later, then we let tk — )• 0"'" so that tt(tfc) ^ ^, and we get 



< 



u{t' ,x) dx 



(7.6) 



Then we put t' = t'^ and let — >• 0"*" so that u{t'^) w(t^) — )• (^(0) = and we obtain the first 

inequality 



99 d/i 



< 



93 du 



(7.7) 



Then, we proceed exactly in the same way but we exchange the roles of and t'^ to obtain the opposite 
inequality [ Jj^d ip d/i] ° < [ j^d ^ dv\ " . Therefore we have that [i = v bs, positive linear functionals on 
(M'^) as desired. □ 

Theorem 7.2 (Existence and uniqueness of initial trace, FD case) Let < m < 1 and let u 

he a nonnegative weak solution of equation (1.1) in {0,T] x M"^. Assume that \\u{T)\\i^i(^^d-^ < 00. Then 
there exists a unique nonnegative Radon measure fi as initial trace, that is 



ijj d/i 



lim 



n(t, x)'4}{x) dx , for all ip G Co( 



Moreover, the initial trace ^ satisfies the bound 

li{BR{x,)) < ||n(r)||Li(M.) + 

where Ci = Ci{m, d, s) > as in ( 2.3 ) . 
Proof. The proof is divided into three steps. 



(7.8) 
(7.9) 



• Step 1. Weighted estimates I. Existence. First we recall the weighted estimates of Theorem 2.2 
which imply for all < t < Ti < T 

l—m / f \ 1— m 

ti(t, j;)0/j(x) dx I < 



u{T, x)(I)r{x) dx + |r - Ti I 



(7.10) 



since < 1 and where Ci > depends only on a, m, d as in Theorem 2.2 . Since (fiR > 1 on Br it is 
clear that this implies hypothesis {i) of Lemma 7.1, therefore it guarantees the existence of an initial 
trace that satisfies the bound ^i{Br{xo)) < Ki = ||n(r)||Li(Kd) + Cii?'^(i-"^)-2« T . 

• Step 2. Pseudo-local estimates. Uniqueness. In order to prove uniqueness of the initial trace is is 
sufficient to prove hypothesis {ii) of Lemma 7.1, namely we need to prove that 



n(t, x)iIj{x) dx 



< 



u{t' , x)ip{x) dx 



+ K2L0{\t-t'\) 



(7.11) 



for aU < t,t' < Ti < T and for all t/j G C^{R'^) . We will see that this is true for a = 1 and 
uj{\t - t'\) = \t- t'\. Let V G C^iW^), then we have 



d 
dt 



u(t, x)ip{x) dx 



< 



(-A)*u'"Vdx 

1 

\{-Amx)\ 



--{a) 



< h \\<t)R\\hHm Ki 



-^Yipdx 

1— m 



< 



(j)R dx 



U(j)R dx 



K' 



2 • 
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Notice that in (a) we have used the fact that (—A)* is a symmetric operator. In (b) we have chosen 



(j)R{x) := 4>{x/R) , with (p as in (2.2) of Lemma 2.1, with the decay at infinity a = d + 2s. It then 
follows that 



|(-A)>(x)| 



<k7, 



since we know by Lemma 2.1 that |(— A)*^/^(x)| < ^^^-(d+2s) ^ have chosen (j)p> > kQ/\x\'^~^'^'^ . 

We have also used the fact that the L^-norm (m < 1) is less than the L^ norm since the mea- 
sure (pFtdx is finite. In the last line of the display we have used the bound of Step 1, namely that 
{J^au{t,x)ipR{x)dxy~"' < Mt + Cii?'^(i-")-2s^ ._ ah < t < Ti . Summing up, we have 

obtained: 

u{t, x)ip{x) dx 



dt 



< K2. 



Integrating the above differential inequality we obtain: 

u{t,x)'4j{x)dx < / u{s,x)'4j{x)dx + K2\t- s\ ioi any s,t>Q a.nd aW i) ^ C^{W^) . (7.12) 
Step 3. We stih have to pass from test functions V G C^(M'^) to V G C°(M'^) in formula ([7^, but 



this is easy by approximation (mollification). □ 

Remarks. (i) The proof applies with minor modification to the class of solutions with data uq € 
'L} {W^ , dx) constructed in Section 



(ii) Notice that estimates (7.12) are only pseudo-local estimates: the global information about u{T), 
namely the bound ||'u(r) ||Li(iRti) is contained in the constant Ki and therefore in K2. 

(iii) The existence of solutions and traces for the standard FDE with (not necessarily locally finite) 
Borel measures as data is studied in Chasseigne- Vazquez [7]. We do not address the corresponding 
question here. 

Theorem 7.3 (Existence and uni queness of initial trace, PME case) Let m > 1 and let u be 

a solution to the Cauchy problem 1.1 on (0,T] x M^. Assume that \\u{T)\\j^i^^d'j + \\u{T)\\i^oo^^d'^ < +oo. 



Then there exists a unique nonnegative Borel measure /i as initial trace, that is 

ipdfi= lim / u{t, x)'ip{x) dx , for all ip £ Co(R'^) 
Moreover the initial trace fi satisfies the bound 

( j^2s+d{m-l) ' 



KBr{xo)) < Ci 



1 

m — 1 



+ T'is u{xo,T)2si 



(7.13) 



(7.14) 



where Ci = Ci{m,d,s) > as in Theorem 6.1 



Proof. The proof is divided in three steps 



in the form (6.3 ) 



Step 1. Weighted estimates L Existence. First we recall the lower bounds of Theorem (|6.1|) rewritten 

' j^2s+d{m-l) " 



Br{xo) 



u{s, x) dx < Ci 



1 

m — 1 



+ T2s n(xo,T)2sti 



(7.15) 
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on the time interval (s,T] C (0,T] . It is clear that this implies hypothesis (i) of Lemma 7.1, therefore 
it guarantees the existence of an initial trace that satisfy the bound ii{Bii{xq)) < Ki . 

• Step 2. Smoothing effects and mass conservation. In |12j there are the global L-*^ — smoothing 



effects which provide global upper bounds for solutions to the Cauchy problem 1.1 . We apply such 



smoothing effects to solutions to our reduced Problem 4.5 to get 



I'"(*)IIl°°(R'«) 



< 



^ Jo 

1^ 



(7.16) 



where = l/[2s + d{m — 1)] and the constant /qo only depends on d,s,m. Moreover, we know that 



there holds also the conservation of mass on the time interval [t/2, T] C (0, T], so that inequality (7.16) 
becomes 

||n(t)||Lo.(K.)<^||n(r)||2l'J^,). (7.17) 



Step 3. Weighted estimates II. Pseudo-local estimates. Uniqueness. In order to prove uniqueness of 



the initial trace is is sufficient to prove hypothesis (ii) of Lemma 7.1, namely we need to prove 



n(t, x)ip{x) dx 



< 



u{t' , x)ip{x) dx 



+ K2L0{\t-t'\) 



(7.18) 



for ah < < Ti < r and for all -0 G C^{R'^) . We will see that this is true for a = 1 and 
a;(|t -t'\) = It" - t'^l) with a = 2s/[2s + d{m - 1)] . Let e C^{R'^), then we have 



d 
dt 



u(t, x)tl'{x) dx 



-A)"M'"Vdx 



= (a) 



u'^i-Ayijdx 

u{t) dx 



< u 



')d-d(m—l) Tin— I 



rfd-Sim—l) rm— 1 



IL°° 



IL°° 



I^(*)IIl°°(R<') 
•= 

■j-d'd{Tn—l) 



Notice that in (a) we have used the fact that (—A)** is a symmetric operator. In (6) we have used the 
smoothing effect ( |7.17 ) of Step 2 . Summing up we have obtained: 



dt 



u{t, x)ip{x) dx 



< 



K2 



ld-d{m-l) ■ 



Integrating the above differential inequality we obtain for any s, t > 0: 



/ u{t,x)4>{x) dx < / ^(5, x)'0(a;) 



dx + 2519 K2 



for all G C~(M'^) . (7.19) 



• Step 3. We stiU have to pass from test functions V G C^(M'^) to G C°(M°') in formula ( |7.13 ) 
but this is easy by approximation (mollification). □ 



We notice that the estimates (7.19) are only pseudo-local estimates: the global information about 



u{T), namely the bound ||'u(r) ULimd) is contained in the constant K2. 
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8 Appendix: Complements and computations 
8.1 Definition of the fractional Laplacian. 

According to Stein, [15], chapter V, the definition of the nonlocal operator (— A)*^/^, known as the 
Laplacian of order a, is done by means of Fourier series 

((-Ar/Vr(x) = (2vr|x|r/(x), (8.1) 

and can be used for positive and negative values of cr. If < a < 2, we can also use the representation 
by means of an hypersingular kernel, 

{-Ar/'g{x) = c,,.P.V. [ f dz, (8.2) 

jRd \X- 

where Cd^a = ^ Tr'^/^^\^i~^J/2^^ ^ normalization constant. Another classical way of defining the fractional 
powers of a linear self-adjoint nonnegative operator, in terms of the associated semigroup, which in our 
case reads 

1 , ,A , . , dt 



{-Ar/'g{x) = — — / (e*^5(^) - 9{x)) (8.3) 

J- (-2)^0 

In this paper we consistently put a = 2s, < s < 1 (sometimes, also s = 1). 
8.2 Definition of weak and very weak solutions 

We recall here the definitions of weak and strong solutions taken from [T^]. We finally introduce the 
definition of very weak solutions. 



Definition 8.1 A function u is a weak solution to Equation (1.1) if: 

• u£ C((0,oo) : Li(M^)), |n|'"-iu € L2^^((0,oo) : H%R'^)); 

• The identity 



Jo Jr'' ot Jq j^d 

holds for every ip G Cg(M"' x (0, 00)); 



u^dxds-T f {- Ay /Wur-^u){- Ay/^if dxds = 0. (8.4) 



A weak solution to Problem (1.1)-(1.2) is a weak solution to Equation (1.1) such that 



u G C([0, 00) : Li(M'^)) and u(0, ■) = uq e L\R'^) 



moreover 



Note that in |12j these weak solutions are given the more precise name weai: -energy solutions. We 
recall that the fractional Sobolev space if*(M'^) is defined as the completion of C^(M'^) with the norm 
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Definition 8.2 We say that a weak solution u to Problem (1.1)-(1.2) is a strong solution if moreover 
dtu G L°°((t,oo) : L^{W^)), for every r > 0. 



Definition 8.3 A function u is a very weak solution to Equation (1.1) if: 



u e C((0, oo) : L\^^{W^)), \u\^-'^u e L\^^ ((0,oo) : (M^ (1 + dx)); 
The identity 

I I u^dxds-[ [ |np-in(-A)V dxds = 0. 
holds for every G C^([0, T] x M'^) ; 



(8.5) 



A very weak solution to Problem (l.l)-0.2^ is very weak solution to Equation (1.1) such that 



moreover u £ C([0, oo) : Lj^^iM.'^)) and u{0, ■) = uo € Uod^'^)- 



8.3 Some functional inequalities related to the fractional Laplacian 

We recall here some useful functional inequalities which have been used throughout the paper. 

Lemma 8.4 (Stroock-Varopoulos' inequality) Let 0<s<l,g>l. Then 

4{q-l) 



\v\'^''^v{-Ayv dx > 
for all V G L9(M"') such that {-AYv G L'2(M'^) 



(-A)2|?;|2 



(8.6) 



Remark. We have used the above Stroock-Varopoulos inequality, applied to < v = and q 
{p + m — 1) /m > 1 , whenever p > 1, which is 



u\P-\{-Ay{\ur^'u) dx > 



m— 1„ 



Am{p — 1) 



{p + m — 1)^ 

Theorem 8.5 (Sobolev Inequality) Let < s < 1 and 2s < d. Then 

11/11^ < 5, (-A)^/2/ 

d-2s 

where the best constant is given by 



, s p+m — 1 
A)2|ti| 2 



dx . 



(8.7) 



C2 . r) — 2s „— s 

o„ := z vr 



r(^) 



T{d/2) 



, 2. 



(8.9) 



and is attained on the family of functions 
F{x) := a [6^ + (x — xq 



d — 2s 

21 — 



with X, Xq G M"" and a G M , b> . 
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8.4 Proof of Lemma 12.11 



Proof. The proof is divided into several steps. 

• Step 1. The integral is convergent. First we have to prove that 



cT1|(-A)V(x) 



ip{x) - Lp{y) 



\x-y 



d+2s 



dy 



< oo for any x G 



to this end we fix x G M'^ and we split the integral in two parts: 

/ 



V{x) - y(y) . 

\d+2s 



F-yr^ J\x-y\>5 \X-y\'^^-' J\x-y\<5 \X - y 

where 5 > is taken so small that the following Taylor expansion around x € holds true 

ip{y) = ip{x) + Vip{x) ■{y-x) + {y - x)* T)'^ip{x) {y - x) 
for some x ^ Bi (x) . Therefore we have 

ip{x) - (f{y) 



\d+2s 



^ dy 

\x-y\<S \X ~ y 

V(^(x) • (y - x) 



\x-y\<5 \X — y 



d+2s 



d2/ + 



(y - xyV)^Lp{x) {y - x) 



\x-y\<& 



<{a) sup ||9jjV3||L^(Kd) 
l<i,j<d 



l\x 



l<i,j<d 

where in (a) we have used that 



\x-y\<s \x-y 
dr 



d-(2-2s) 



x-y 
dy 



d+2s 



dy 



(52(1-^) 

< sup ||a,jy.||Loo(iRd) ^i_2(i-s) =W ^(2(1 



P.V. 



V^ix) • (y - X) 

d+2s " 



\x-y\<5 \X — y 



for symmetry reasons. In (b) we used the fact that \dij(p{z)\ < K for some positive constant K that 
depends only on a . On the other hand, the outer integral is easily seen to be finite, indeed 



\x-y\>5 \x y 



(p{x) - ip{y) 
dT2rdy 



< 2||v5||l° 



1 



\x-y\>5 \X y 



d+2i 



dy 



< 2ujd 



dr uJd 



l+2s g^2s ■ 



The above estimates for / and // do not depend on x G M*^, hence |(— A)'*(^(x)| is finite for all x G M'^ . 

• Step 2. Better estimates for \x\ large. We are going to use the hypothesis that (p is radially 
symmetric and decreasing for |x| > 1 and that (p{x) < |x|~" , \D'^ip{x)\ < co|x|~"~2 , for some positive 
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Figure 3: The 4 regions in which we split the integral 



constant a and for |x| large enough. We are interested in the behaviour of |(— A)*(/9(2;)| for large values 
of X, therefore we fix x G M*^ with |x| sufficiently large. We have to estimate 



c7i|(-A)V(x)| 



ip{x) - ip{y) 



to this end we split the integral into four parts, see Figure [3| 

ip{x) - tf{y) 



ip{x) - ip{y) 

1 \d^'^y 

\x — 



+ 



d+2s 

d+2s dy + 



\y\>3\x\/2 \X-y 
<f{x) - (p{y) 



dy + 



{|x|<2|s,|<3|x|}\B|,|/2(x) \X-y 
(fix) - <f{y) 



<f{x) - (p{y) 

dT2rdy 



^.10) 



\y\<\x\/2 \x-y 



\d+2s 



dy = 1 + 11 + 111 + IV 



We estimate the four integrals separately, keeping in mind that we are assuming (/3 > in this latter 
case. The first integral can be estimated as follows 



I 



\y\>3\x\/2 \X-y 



dr 

1 



k 



\a+2s 



since ip{y) < ip{x) when \y\ > 3\x\/2, therefore \(p{x) — (p{y)\ < f{x) , and we remark that the constant 
ki depends only on a,s,d, since ip{x) < \x\~" and \x\ is large enough. The second integral gives 



// < 



(p{x) - ip{y) 



{\x\<2\y\<Z\x\}\B^^u^{x) \X-y 



d+2s 



dy 



< 



^ix/2) 



{\x\/2) 



d+2s 



3N 
2 



kn 



\x 



a+2s 



since </?(y) < ip{x/2) when \y\ > |x|/2, therefore \^{x) — '^{y)\ < ip{x/2) , and we remark that the 
constant /c2 depends only on a, s, d, since ip{x/2) < \x\~" and |x| is large enough. 



33 



We can estimate the third integral as follows: 

Lp{x) - ip{y) 



III 



\d+2s 



dy 



Vv?(x) ■ {y- x) 



\d+2s 



dy + 



{y — x)* Hess(/3(x) (y 



<(a) sup ||9ij,-(/'||L°°(B|^| (x)) 
l<ij'<d 



|d-(2-2s) 



x-y 
dy 



d+2s 



dy 



< 



\ct+2 



dr 



j,l-2(l-s) 



|«+2 



2-2s 



a+2s 



where in (a) we have used that 



P.V. 



Vf{x) ■ (y-x) 



\x~y\<5 



\d+2s 



dy = 



for symmetry reasons as in Step 2. In (6) we used the fact that — x| < \x\/2 implies \x\/2 < \z\ < 
3|x|/2, therefore \dfjip{z)\ < co/\z\°'~^'^ < 2°+^co/|x|"+^ for all z e B\xy2{x), recalling that |x| is always 
taken large enough. The constants k'^ and depend only on a, s, d. 
It only remains to estimate the fourth integral: 



IV < 



\ip{x) - ip{y)\ 

\d+2s 



2d+ 



2s 



\X 



d+2s 



'\y\<\x\/2 \X 

since we observe that \y\ < \x\/2 implies < (p{2y) < (p{y) which gives \(p{x) — if{y)\ < fiy), 

moreover we also have that |y| < |x|/2 implies that |y — x| > |x|/2. The term represents the long-range 
influence of the inner core of the function at large distances and will make for different conclusions of 
the lemma depending on the case. Indeed, we have the following estimates for \x\ large enough: 

• If a > d the last integral is finite and we get IV < k^/lx]'^'^'^^ . 

• If a < d the last integral grows like Ixl*^"" and we get IV < k^/\x\°'~^'^^ . 

• Finally when a = d we get IV < fcg log |x|/|x|'^"'"^* . 

We finally remark that the constants k4, k^, kg depend only on a,s,d . 

• Step 3. Positivity estimates for |x| large. In the case when a > d we need to prove that if (/? > 
then we have that |(— A)^99(x)| > C4|x|~('^"'"'^''^ for all |x| > |xo| >> 1 . We split the integral into four 
parts, as in Step 2, equation (8.10), see Figure [sj 

fix) - f{y) 



-^d,s 



^(-A)V(x) 



Rd \x-y 



d+2s 



dy = 1 + 11 + 111 + IV 



We have proven in Step 2 that |/| + + < (/^i + ^2 + /cs) / Ixl""*"^^ and we recall that the constant 
ki depend only on a, s, d. We just have to obtain better estimates for the last term, to this end we 
further split the integral in two parts: 



IV 



ip{x) - ip{y) 



\y\<\x\/2 \x-y 



d+2s 



dy 



ip{x) 



\y\<\x\/2 \x-y 



d+2s 



dy 



\y\<\x\/2 \x-y 



d+2s 



dy = IVa - IVb 
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Let us calculate 



0</F„ 



KN/2 



d+2s - 



k - y| 



{\x\/2) 



d+2s 



KN/2 



A;4 



\a+2s 



since |x — y| > |x|/2 when \y\ < \x\/2 and ip{x) < \x\ We remark that the constant depends only 
on a, s, d. On the other hand, /VJ, > and 



< IVb 



dT2ldy < 



1 



j<N/2 

Summing up, we have obtained that 

ip{x) - <f{y) 



\<\x\/2 



\d+2s 



-(-A)V(x) 



\x - y\ 



d+2s 



dy > Cd,s [IVb + \II\ + + \IVa\ )] 



\d+2s 



\a+2s 



/C5 



Cd,s ^ C4 



since |/| + \II\ + |///| + \IVa\ < (ki + fcs + fcg + k^) = /.^ and C4 > since a > d, if we 

choose \x\ sufficiently large, namely |x|"~'^ > fe5/||93||Li(-]gd) . □ 



8.5 Optimization Lemma 



We state and prove here a simple technical lemma that has been used in the proof of Theorem 4.1 



Lemma 8.6 Let < m, s < 1 , 2s > d{l - m) , 1!) = l/[2s - d{l - m)] > and B,C,t > 0. Define 

A{t) Bt _ C 



F{t,R) :-- 



J^d{l-m) Jl2s 



with A{t) := M 



fd{l-m)-d ■ 



Then there exists 



t,:=2s^{—j >0 



and 



R{t) 



2sBt 



J{1 - m)A{t) 
so that for all t>t^, we have 

F(R{t),t) = 



2s 

d(l — m) 



2s 



(2si?) 



d{l - m) (2s??)<^(i-™) - 1 



j\^d{l-m) 



(8.11) 



> (8.12) 



d{l — m) 
27 



-I 2s'd 



2sS 



> 0. 



Proof. First we observe that A(t) is monotone increasing in t > 0, and that ^(t*/2si?) = 0, where f^, 



has the expression given by (8.11), so that A{t) > ^(t*) > A{t/2s'd) = since 2s'!? > 1, and 



A{t^) = ^ , -^M > 
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Moreover, it easy to check that t^, is also the value for which A{t^) — t^A'{t^:) = 0. Next we fix a time 
t > t=K and we find the maximum with respect to R of the function F{t, R): 

a Flfl f) = dil-m)A{t) 2sBt 
and diiF{R{t),t) = 0, so that the maximum is attained at R{t) whose expression it is easily checked 



to be (8.12). It only remains to prove that R{t) > R{t*) > 0, to this end we observe that 

' ^ Y'^ A{t)-tA'{t) 



dtR{t) = I? 



2sB 



d{l — in) 



t 



Ait) 



Ait) 



2 



and it is clear now that the minimum is attained at , since (9ti?(t*) = 0, because we already know 
that Ait^) -t^A'it^) = 0. □ 

8.6 Reminder about measure theory 

We recall here some basic facts on measure theory for convenience of the reader. We refer the interested 
reader to the books [9l [T^. 

Definition 8.7 A measure /U is regular if 

\/A C M"^ 3B ^-measurable such that A Q B and niA) = fJ-iB) . 

A m,easure fx is Borel if every Borel set BiW^) is ^-measurable. A measure fi is Borel regular if 

VA CR'^BB e BiR"^) such that A C B and fiiA) = fiiB) . 

A measure ^ is Radon if is Borel regular and niK) < +oo for any compact set K C . 
A sequence of measures Hn converges weakly (star) to the measure ji, ^ l^ as n ^ oo if 

lim / ipdi^n= 9?d/i for all cp G C°(M'^) . 



n— >oo 



Theorem 8.8 (Weak compactness for measures) Let {fJ-n} be a sequence of Radon measures on 
R'^ satisfying 

sup UniK) < OO for any compact set K CR'^. 

n 

Then there exists a subsequence fin^ and a Radon measure /i such that ^ as A; oo . 

Theorem 8.9 (Riesz Representation Theorem) Assume L : C^(M'^) — )• M is linear and nonneg- 
ative, so that 

Lip>0 for allO<ipe Cf 
Then there is a unique Radon measure ji on R'^ such that 

Lip = ipd^ for all (p E 
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